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1. Introduction 

During last decade, the attention of many authors was attracted to quantum phenomena 
in cavities and media with moving boundaries, known under the names non- stationary 
Casimir effect (NSCE) dynamical Casimir effect 0, or mirror (motion) induced 
radiation For the most recent achievements in this field and references to 

other works see, e.g. [5-23] (the problems of the classical electrodynamics with moving 
boundaries were studied, e.g. in |2^, |^). A complete analytical solution to the problem 
of a one-dimensional (Fabry-Perot) ideal cavity with resonantly vibrating boundaries 
was found recently in [0. This solution holds for any moment of time (provided the 
amplitude of the wall vibrations is small enough); moreover, it enables not only to 
calculate the number of photons created from an arbitrary initial state, but also to 
account for the effects of detuning from a strict resonance. 



In the present article, continuing the line of research of ||T^, we calculate the effects 
of squeezing and find the photon distribution function of the field in each mode inside 
the cavity. Although a possibility of squeezing the electromagnetic field in a cavity with 



moving boundaries due to the NSCE was discussed for the first time ten years ago ||2^ 



the concrete calculations of the variances of the field quadrature components were made 



only in the short-time [^, and long-time ^ limits under the condition of the 
strict resonance. As to the photon statistics in an oscillating Fabry-Perot cavity, it was 
not discussed at all until now (the photon distribution in a 3D nondegenerate resonantly 
driven cavity was found in 0). 

2. Field operator in a ID cavity with oscillating boundaries 



We start with a brief description of the main results of |17| which are used in the 
subsequent sections. Following the model of a "scalar electrodynamics" |3^ , we assume 
that a ID cavity is formed by two infinite ideal plates whose positions are given by 
xieft = and 

Xright = L{t) = Lo {I + e sin [puji{l + 6)t]) , |e| < 1, |5| < 1 (2.1) 
where Ui = ttc/Lq and p = 1,2, . . .. As was shown in |jT^, the results can be easily 



generalized to a generic case when both the mirrors oscillate (see also |11, ^). We are 
looking for the only component of the operator vector potential of the electromagnetic 
field A{x,t) in the Heisenberg representation in the form (hereafter c = h = l) 

oo n 

A{x,t) = E ^ \k'^^''\x,t) + h.c.l , \k,bl] = 5nk. (2.2) 



n=l 



This operator must satisfy the wave equation Att — A^x = and the boundary conditions 
||30| A(0, t) = A{L{t),t) = 0. For t < 0, when the wall is assumed to be at rest, function 
tlj^"'\x,t) has a simple factorized form il)Q^\x,t) = e"*"^"* sin (vrnx/Lo), oOn = nuoi. The 
normalization factors 2 / -s/n in (|2.2|) are chosen in such a way that the energy of the 
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field in the stationary case can be represented as a sum of energies of independent mode 
oscillators (see the next section). 

For t > 0, following |^ we expand each function ^^^\x,t) in a series with respect 
to the instantaneous basis 

irkx 



^^''\x,t) 



(2.3) 



This way we automatically satisfy the boundary conditions. Then the wave equation is 
replaced by the infinite set of coupled differential equations 



(2.4) 



where 



UJkit) = kTc/L{t), gi^j = -gj^ = (-1) 



k-j _ 



2kjLit) 



{f-k^)L{ty 

As was shown in the set of equations ( p.4| ) can be significantly simplified in the 
case of the (quasi)resonance oscillations of the wall given by the dependence (|2.1|) , if 
one writes 

g(")(t) = pWe— ^(1+'^)* _ pWe-^(i+'5)* (2.5) 

assuming the coefficients p^^^ {k = ±1,±2, n = 1,2,...) to be slowly varying 
functions of time, whose derivatives are proportional to the small parameter e. Putting 



]5D into equation (|2.4| ) and neglecting the second order terms like p\. 



one 



obtains, after averaging over fast oscillations with the multiple frequencies of uii, the 
equations [0 



d („) 

-rPk 
dr " 



a 



{k + v)ptlp -{k- v)pTlp + 2i7^P, 



(n) 



(2.6) 



where 7 = 5/e, a = (—1)^, and r = ^euit is the 'slow time'. The exact solutions to 



equation (|2.4| ) differ from the approximate form (|2.5|) by the terms proportional to the 
higher harmonics exp [±iruJkt], r = 2, 3, . . .. However, the magnitudes of the corrections 
are of the order of e (or less), so they can be neglected, at least under the condition 
t <^ ^2 ~ (^1^^) • Foi' the values Ui ~ 10^° s~^ and e ~ 10~^ corresponding to the 
possible experimental realisations |^ the limiting time t2 is of the order of weeks. 

Due to the initial conditions p^"''' (0) = 5kn the solutions to ( |2.6| ) satisfy the relation 
Pj+mp^ = if j 7^ /c. The non-zero coefficients p!^^ read 

r (1 + n + j/p) (aK)"-™A'"+"+2j/P 



(J+np)( ^ 
rj +mp \ I 



r(l + m + j7p)r(l + n 



X F (^72 + Hp , — m — j/p ; 1 + n — m; K 

where F(a, b; c; z) is the Gauss hypergeo metric function, 

sinh (apr) 



+ sinh^ (apr) 



7 



A 



[2.7) 



(2.8) 
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The functions ( p.7| ) are exact solutions to the set of equations ( |2.6| ) relating the 
coefficients with different lower indices. Besides, these functions satisfy another set 
of equations, which can be treated as recurrence relations with respect to the upper 
indices 

d 



dr 
d 

dr 



n < a 



{■- 



Jn-p) _ Jn+p) 
rm rm 



' (p-n)* _ (p+n) 
r—m rrri 



+ 2^7P^::)}, n>p, p(^)^0 



(2.9) 



+ 2z7p(;:)}, n = l,2,...,p-l (2.10) 



The consequences of equations (p.6|), (|2.9|) and (|2.10|) are the identities 



m=—oo 



o 

E 

1 r 



nk 1 ri, k 1, 2, 



n=l 



^ yrm rj r—in H—j 



■'mj ; 



m,j = 1,2, 



oo 1 



0, m,j = l,2. 



(2.11) 
(2.12) 
(2.13) 



The formulas given above hold for any value of the detuning parameter 7. For 
7 > 1 one should replace the functions sinh(ax)/a and cosh(ax) by their trigonometrical 
counterparts sin(ax)/a and cos(ax), where d = " 1- 



3. Squeezing the initial vacuum 

We suppose that after some interval of time T the wall comes back to its initial position 
Lq. For t > T, the field operator assumes the form 



A(x,t) = y ^ sin (vrnx/Lo) fa„e-^'^"(*+^^) + h.c. 
^1 ^ 



(3.1) 



where operators related to the initial operators 6„ and b}^ by means of the 

Bogoliubov transformation (r^ = ^euiT) 

jm 
n 



E [bnpt' M - Pjt {tt) 

n=l ' " 



m = 1,2 



(3.2) 



The commutation relations 



dn ; dj^ 



6nk hold due to the identities ( p.ll| )-( ^T^ ) which 
are nothing but the unitarity conditions of the transformation (|3.2|) . These commutation 
relations together with the expression for the energy of the field 



1 fLo 

H = — dx 
Svr Jo 



dt ) ^[dxj 



CO , 

Y.^n[didn^ - ] (3.3) 



convince us that namely dn and d\ are true photon annihilation and creation operators 
at t > T (like the operators 6„ and 6jj were 'physical' ones at t < 0). 
Now we introduce the Hermitian quadrature component operators 

= (dm + d^^ /\/2, Pm= (dm - / (« V^) 
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and calculate their variances Um = ~ i'lm)'^, Kn = {Pm) ~ {Pm)'^ and the covariance 
= ^{PmQm+fjmPm) — {Pm) {(jm) in the vacuum State with respect to the initial operators 
bn- bn\0) = (remember that we use the Heisenberg picture). Using (3^) we obtain 
(assuming hereafter ui = 1) 

U =!!^y ' \/ =-y + ' (34) 



2 fr^ n 



n=l " n=l 



Y^ = Y: -Im [pL"^V^1J (3.5) 

n=l 

where the coefficients p|^2i should be taken at the moment T, thus their argument is r^. 
Strictly speaking, the expressions (|3^)- (|3l5| ) have physical meanings at those moments 
of time T when the wall returns to its initial position, i.e. for T = A^7r/[p(l + 5)] with an 
integer N . Consequently, the argument tt of the coefficients in (p.4|)-( p75D assumes 



discrete values r*^^-* = N ei: / \2p{l -\- 5)]. One should remember, however, that something 
interesting in our problem happens for the values r ~ 1 (or larger). Then ~ ^ 1, 
and the minimal increment Ar ~ e is so small that tt can be considered as a continuous 
variable (under the realistic conditions, e < 10^® [^]). For this reason, we omit hereafter 
the subscript T, writing simply r instead of Tt or t^^\ 

Differentiating the right-hand sides of equations ( |3.4| ) and ( |3.5| ) with respect to the 
'slow time' r, one can remove the fraction 1/n with the aid of the recurrence relations 
(2.9) and ( p. 10 ). After that, changing if necessary the summation index n to n±p, one 
can verify that almost all terms in the right-hand sides are cancelled, and the infinite 
series are reduced to the finite sums: 



dK^/dr 



p-1 



am 

n=l 

p-1 



(p-n) (v-n) 
rm h'—m 



h'—m Km 



dFm/dr = am E Im {pt^'' P^-""^* + P^-lp^^m^) 

n=l 

Now one should take into account the structure of the coefficients ( |2.7| ): they are 
different from zero provided the difference between the upper index n and the lower 
one m is some multiple of the number p. If m = j ' + pk with j = 1, ... ,p — 1 and 
/c = 0, 1, 2, . . ., then only the terms with n = j or n = p — j survive in the sums above. 
Depending on whether j = p/2 or j ^ p/2, we obtain two different sets of explicit 
expressions for the derivatives of the (co)variances. 

1) If m = J ' + pk but j 7^ p/2 (in particular, for all odd values of p), then 

§^ = f^ = 2.mRe(pijVLV'). ^ = 0. (3.6) 

In this case = and Um = Kn = A/'m + 1/2, where Afm is the mean number of photons 
created from vacuum in the mth mode calculated in |T^, so there is no squeezing. For 
7 < 1 the quadrature variances monotonously increase in time, with an asymptotical 
linear dependence. If 7 > 1, the variances oscillate in time with amplitudes inversly 
proportional to 7^ — 1, being always not less than 1/2. 
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2) Some squeezing can be achieved only in the 'principal' modes with the numbers 
^l = p{k + 1/2), A; = 0,1,2,...: 



— /LiRe 



dr 

dy^/dr = /ilm 



(p/2) _ (p/2) 



dr 



/iRe 



+ 



{p/2) 



P-^ 



(3.7) 
(3.8) 



(in particular, it is necessary that p be an even number). In the strict resonance case 
(7 = 0) all the coefficients p^^^"^^ are real, so 1^ = and dU^/dr < in the whole interval 
< r < 00, resulting in the inequality Uf^i^r) < 1/2. 

Note that the coefficients Pp^+p/2^ depend on the parameter p only through the 
dependence of the the variable n on the product pr: see equation ( p.8| ) and the explicit 
form of these coefficients in Appendix A . Thus it is sufficient to consider the most 
important special case of the parametric resonance at the double fundamental frequency 
2uji (i.e. p = 2), since the formulae for p > 2 can be obtained by a simple rescaling 
of the 'slow time' (for the 'principal' modes). In this case, only the odd modes can be 
excited from the vacuum, and they do exhibit some squeezing. 

Using equations ( |3^ ) and ( p.8| ) one can immediately find the Taylor expansions of 
the (co)variances at r — (assuming (—1)!! = 1): 




2m+l 



(2m 



m! 



It 



2m + 1 
(m + 1)2 



2m+l 



-27(2m+ l)r 



2(m+l) 



(2m 



m! 



+ 



(3.9) 



(3.10) 



We see that the [/-variances are always less than 1/2 at the initial stage, but the degree 
of their squeezing rapidly decreases with increase of the number m. Note that the 
dependence on the detuning parameter 7 in the short-time limit appears only in terms 
of the order of r2'^+^ (and higher). 

In the opposite limit r — 00 (or k — 1), using equations (|3.7] ), (|3.8| ) and the 
asymptotics of the Bogoliubov coefficients (|A.4| ) we obtain constant time derivatives 



d[/'2m+i/dr|^ 
dV2m+i/dr|^_ 
dl^2m+i/dr|^_ 



16a 



7r2(2m + 1) 

16a 
7r2(2m+ 1) 
8a 



sm 



cos 



1 

m H — 
2 

1 

m H — 
2 



(3.11) 
(3.12) 

7r2(2m+l)^"^^^-^ ' ^^"^ ^^-^^^ 

where (j) = arcsin7. Consequently, all the (co)variances increase with time linearly, 
giving the constant photon generation rate in the 'principal' (odd) modes 

3a 



■ sin [(2m -|- 1) 1 



dA/'2m+l/d7 



(3.14) 



7r2(2m+l) 

in agreement with fl^. Equation ( p. 141) results in a simple estimation of the mean 
photon number in the /ith mode at r > 1: A/'^(r) ^ ar/fi. 
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Since the covariance is different from zero if 7 7^ 0, tlie initial vacuum state 
of tfie field is transformed to the correlated quantum state One should 

remember, however, that the values of and 1^ yield the (co)variances of the 

field quadratures only at the moment t = T (when the wall stopped to oscillate). At 
the subsequent moments of time the quadrature variances exhibit fast oscillations with 
twice the frequency of the mode. For example (omitting the mode index), 

ag{t') = Ucos'^iut') + Vsin^{ujt') + F sin(2cut'), t' = t-T. 

Therefore the physical meanings have not the values f/^, V"^ and themselves, but 

'a values of the quadrature variances 



rather the minimal 



Ufj, and maximal amax 



during the period of fast oscillations 



(3.15) 



and 



Only in the special case of the strict resonane (7 = 0) we have 

The equations ( p.7|) and ( p.8|) can be integrated for any values of and 7 in terms of 
the complete elliptic integrals K(k) and E(k): see [Appendix B| for the technical details. 
In particular, for the fundamental mode we find 

^ ' 2(/5-/€)K(/€)E(fi:)+/3E2' 



2{I3 + K)K(K)E(/t) - k^{l3 + k)'K^{k) - /3E2(k) 



27 
7r2 



2K(k)E(k) + E'(k) 



where (3 = ReA = \/l — 'y'^K-^ and k = — k?- Then equation ( p.l5|) yields the minimal 
and maximal invariant variances 
2 



Ml 



2(1 - k)K{k)E{k) + E\k) 



(3.16) 



^^1 



2(1 + k)K(k)E(k) 



(3.17) 



which depend on the detuning parameter 7 only implicitly, through the dependence on 
7 of the function /t(r) ( p.8|) . In the short time limit r <C 1 (then k 2t) we obtain, 
using the Taylor expansions of the complete elliptic integrals, Ui = ^ — r + r"^ -\ — ■ and 



i + r + r2 + . . 



in accordance with [p6l. More precisely, 



2 V 2 4 32 

The minimal variance Ui monotonously decreases from the value 1/2 at t = to the 
constant asymptotical value 2/71^ at r 3> 1, confirming qualitatively the evaluations of 
P^ , and giving almost 50% squeezing in the initial vacuum state. The variance of 
the conjugate quadrature monotonously increases, and for r ^ 1 it becomes practically 
linear function of time: Vi{t ^ 1) ^ IGt/tt^. The asymptotical minimal value 
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Ui{t = oo) does not depend on 7 provided 7 < 1 (only the rate of reaching this 
asymptotical value decreases with 7 as y/l — 7^). In the strongly detuned case, 7 > 1, 
the minimal variance oscillates as a function of r (being always greater than 2/71^), since 
in this case the function k{t) oscillates between —7"^ and 7"^. 

The minimal variance does not go to zero when r 00 due to the strong intermode 
interaction, which results in a high degree of quantum mixing for each mode. Since the 
state originating from the initial vacuum state belongs to the class of Gaussian states 
(see section |]), the quantum 'purity' Xm = Trp^ of the mth field mode (described 
by means of the density matrix p^) can be expressed in terms of the (co)variances 



as 



-1/2 



-- [4 {UmVm — Y^)] Using equations (|3.11| )- (|3.13| ) one can check that 
for r ^ 1 (see [Appendix C|) . For instance, for m = 1 we have (writing 



simply K and E instead of K(/t) and E(/t)) 



Xi 



E^ - k^K' 



-1/2 



(3.18) 



The initial dependence on k is rather weak: x(/t ^ 1) = 1 — -^k + ■ ■ ■. But when 
K 1, X rapidly goes to zero: x(k ^ 1) ~ {8/1^'^) [In (4/k)]^^^^, with dx/dn —00. 

It is worth to mention that in the resonance case there is no effective interaction 
between different modes of a nondegenerate (three-dimensional) cavity possessing non- 
equidistant eigenmode spectrum. In such a cavity the minimal variance asymptotically 
goes to zero for r ^ 00 |T^. 

For the sake of completeness we bring the formulae (found in |T^) for the mean 
photon number in the fundamental mode and the total energy of photons created in all 
the modes (for the initial vacuum state): 



1 

2 



12a2 



sinh^ (par) . 



(3.19) 
(3.20) 



In [[13 it was derived also the formula for the second derivative of the total mean number 
of photons created from vacuum in all the modes 



Hot 



dr2 



(3.21) 



Now we are able to integrate this equation. Taking into account the condition 
d^ /dr = at r = we obtain very simple expression 

Utot = ^K(«:) [K(/€) - E{k)] . (3.22) 



The expressions for the variances in the modes with numbers > 3 are rather 
involved; their general structure is discussed in Appendix B . Here we give only one 
explicit example - the variance f/3 for 7 = 0: 

2 



+ (1- 



kHi 



k) (4 + 10k + 9k/ 



K 



4k^ - Uk^ - 20k 



8) K{k)E{k) 



+ (Ak'^ + 6k^ 



K 



6k + 4 E 



(3.23) 
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The Taylor expansion of the right-hand side of ( |3.23| ) coincides with the expansion ( |3.9| ) . 
The asymptotical value at r ^ cxo equals U^{k = 1) = 38/(97r^) ^ 0.43. We see that 
the squeezing rapidly disappears with increase of the mode number /i. 

The variance V3 can be obtained from ( |3.23| ) by means of a simple substitution 
K — >■ —K. Therefore the mean number of photons in the third mode is given by 



3k^ - 2) K f2E - k^K) + 2 (1 + 



(3.24) 



371^^2 

It is remarkable that despite all three (co)variances ?7^, and linearly increase with 
time at r 3> 1 in the generic case 7 7^ 0, the minimal variance tends to a constant 
value at r — s> 00: see [Appendix C| . 

concerning 



The results of this section confirm completely the earlier conclusions [26 



the behaviour of the quadrature variances in the short-time limit for p = 2, as well 
as the results of approximate numerical calculations performed in |^ (in the same 
limit) for p = 2 and p = 3. In the long-time limit we see a qualitative agreement 
with earlier approximate asymptotical formulae of p5| related to the behaviour of 
the squeezed quadrature component variance (for 7 = 0), namely, that the squeezing 
effect is strongest for the lowest mode and that it decreases with increase of the mode 
number. However, there are differences in the numerical values of the squeezed variances. 
This can be explained as follows. In ^ only the leading terms of the Bogoliubov 
coefficients (analogs of the coefficients p^-*) were found. But the coefficients p^^{k) 
and p^^Ii{k) become the same in the limit k = 1 (if A = 1) and they do not depend 
on the upper index n: see [Appendix A| . For this reason, although the leading terms 
of the asymptotic expansions enable to calculate correctly the number of photons and 
the unsqueezed variances like V2m+i, these terms are canceled in the expressions for the 
squeezed variances like f/2m+i (if 7 = 0). It was found in |2^, ^ that the difference 
Am = | — t/2m+i(T ^ oo) decrcascs as l/(2m + l) for m ^ 1. The numerical integration 
of equation ([3?7[) yields for the product (2m + 1) the values close to 0.2 for 1 < m < 5 
(according to ([3.23[) , 3Ai = 0.2166 . . .). It seems probable that the limit value of this 
product at m — > 00 equals 2/7r^ = 0.2026..., but we did not succeed to prove this 
conjecture analj^ically. 



4. Influence of initial conditions 



For an arbitrary initial state of the field one can write f/m = U^""^ + AU^, where 
U^"''^^ is given by equation (|3.4[ ); similar expressions can be written for Vm and Ym- The 
corrections due to the nonvacuum initial states are given by 




ReE 



"J 



m 



'in) ^ in) 



Jj) ^ ^(J) 



± 



pin) ^ in) 



rm r—m 



(bib,)- (hi) (b,) 
{bnb,)-{bn){b,) 



(4.1) 



AK. 



Im 



n,3 



m 



rm H—m Hm r~m 



{bih)-{bi){h) 
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+ 



■in) U) 
rm rm 



(") (i) ■ 

P—mP—m 



(4.2) 



where the average values hke (bl^bj) are calculated in the initial state. All the corrections 
disappear in the case of the initial coherent state, = If the initial density 

matrix is diagonal in the Fock basis (as happens e.g. for the Fock or thermal states) 
then (bl^bj) = fn^nj {i^n > 0), all other average values in ( |4.1| ) and ( [4.2| ) being equal to 
zero. In this case the double sums are reduced to the single ones: 



n 



pin) 



mJ2 

n 

2m — Im 



An) 



A\4 



E 



n 



pin) + in) 
rm ' y—m 



pin)* pin) ■ 
rm r—m 



(4.3) 
(4.4) 



We see that the initial fluctuations always increase both the variances Um and Vm (for 
the diagonal density matrix). However, asymptotically at r — oo the corrections are 
bounded for the physical initial states having finite total numbers of photons, because 



the coefficients 



in) ± 



P 



in) 



and Im 



in)* in) 
rm y—m 



do not depend on the summation index 



n in this limit: see equations ([A.6|) and ( |A.7| ). Thus we have in the 'principal' /x-modes 
(for p = 2; 4 



arcsin 7) 

AV^"") 
V 

Ay(oo) 



8Z 



X 



2sin2 (/i0/2) 
2 cos2 (/i0/2) 
— sin (/i0) 



2 = E 



fc=0 



2kTT' 



(4.5) 



These expressions are very similar to ( |3.11| )- p.l3| ). Their consequence is the important 
result that in the limit r ^ 00 the minimal variance does not depend on the initial 
state of the field inside the cavity, provided the initial density matrix was diagonal in 
the Fock basis. The proof is given in [Appendix C| . The correction to the mean number 
of photons tends to the hmit AA/J°°) = 8Z/ (vrV)- 



5. Photon distribution 



Now let us turn to the photon distribution function (PDF) /(n) = {n\pm{t)\'n), where 
|n) is the multimode Fock state, n = (ni, n2, . . .), and Pm(t) is the time-dependent 
density matrix of the mth field mode in the Schrddinger picture. At first glance, there 
is a problem, since all the calculations in the preceding sections were performed in the 
framework of the Heisenberg picture. Fortunately, this problem can be easily resolved 
for a special class of Gaussian initial states (i.e. the states whose density matrices, 
or wave functions, or Wigner functions, are described by some Gaussian exponentials). 
This class includes coherent, squeezed and thermal states; in particular, it includes the 
vacuum state which we are interested in here. 

The solution is based on two key points. The first one is the statement ^ that 
the field evolution in a cavity with moving boundaries can be described not only in 
the Heisenberg picture, but, equivalently, in the framework of the Schrodinger picture, 
with a quadratic multidimensional time-dependent Hamiltonian. The second key point 
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is the fact BO, that the evolution governed by quadratic Hamiltonians transforms 



any Gaussian state to another Gaussian state. 

It remains to take into account that the photon distribution function of any 
Gaussian state is determined completely by the average values of quadratures and by 
their variances |3^, which obviously do not depend on the quantum mechanical 
representation. The explicit formulae in the generic case are rather involved, so we give 
them in [Appendix D| . Here we confine ourselves to the most simple case of the vacuum 



initial state, when all the average values of quadratures are equal to zero. In this case 
the generating function is reduced to [^(2;)]"^/^, i.e. it has the same structure as 

the known generating function of the Legendre polynomials Pn{x). After some algebra 
we obtain the following expression for the photon distribution in the mth field mode: 

2 |(2«„. - l)(2v„ - l)f'^ I 4-i.„;.„ - 1 \ 
1(2«„ + 1)(2.„. + 1)|<"«"'^'^' \^(Anl - 1) (4< - 1) j ' 

It depends only on the invariant minimal and maximal variances Um and Vm- Note that 
the argument of the polynomial in ( ^.1|) is always outside the 'traditional' interval (—1, 1) 
(in particular, this argument is pure imaginary if 2um < 1), being exactly equal to 1 for 
the 'non-principal' modes with Um = Vm = -Afm + |, when formula ( p.l|) transforms to 
the time-dependent Planck's distribution 

Only for the 'principal' /i-modes the spectrum of photons is different from Planck's one 
due to the squeezing effect. The first and second derivatives of the generating function 
( p. ID at z = 1 yield the first two moments of the photon distribution (hereafter we 



suppress subscript m) 



n = ^{u + v-l), an = n'^ - {nf = ^ (2u'^ + 2v^ - 1) (5.2) 
which result in the Mandel parameter 

n- /-I u^ + v^-u-v + 1/2 

Q = an n - 1 = ■ . (5.3) 

U + V — 1 

This parameter appears positive for all values of r, so the photon statistics is super- 
Poissonian, with strong bunching of photons (the pair creation of photons in the NSCE 
was discussed in ||^, ^, |^, |l^). In particular, 

Q2m+i{r ^ 0) ^ [(m + l)(2m - l)!!/m!]' T^"'/{2m + 1), Qi(0) = 1, 

whereas Qm ~ Ki(t) > 1 for r > 1 (if 7 <^ 1). 

The analytical expressions obtained in the paper are illustrated in two figures. In 
figure ^ we show the dependences on the universal variable k ( p.8| ) of the minimal and 
maximal variances ui ( |3.16| ) and vi ( |3.17D together with the purity factor xi ( |3.18| ), the 



mean photon number Afi ( 3.19| ) and the Q-factor ( |5.3| ) for the fundamental mode /i = 1. 



An example of the photon distribution in the principal mode = 1 for 7 = is given 
in figure!^, where the 'cavity' distribution is compared with Planck's one corresponding 
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to the same mean photon number. We see no oscillations typical for the pure squeezed 
states [40-42] (excluding a small 'splash' at n = 2), because the field appears in a 
mixed quantum state (the infiuence of quantum mixing on the oscillations of the PDF 
in generic Gaussian states was studied in |^8|, ^ @)- The asymptotics of the photon 
distribution function ( |5.1| ) in the long-time limit are given in [Appendix D . 



6. Conclusion 



The main results of the paper are as follows. We studied the behavior of the 
electromagnetic field quadrature variances in a one- dimensional cavity with resonantly 
oscillating ideal boundaries in the whole time interval < t < cxd for all field modes and 
for any (small) value of the detuning parameter. We have shown that each field mode 
goes to a mixed quantum state due to the intermode interaction (caused by Doppler's 
effect on the moving mirrors). We found that squeezing can be observed only in the 
'principal' modes with numbers p{k + 1/2), where integer p is close to the ratio of the 
wall vibration frequency to the frequency of the fundamental cavity mode. Analyzing 
the infiuence of the initial nonvacuum state of the field we discovered that the initial 
thermal fiuctuations do not affect the minimal value of the quadrature variance (which is 
less than 1/2) in the long-time limit. This result is important from the practical point of 
view, since it shows that certain significant features of the Nonstationary Casimir Effect 



are not sensitive to the temperature (see also |2^ in the case of the three-dimensional 
cavity). We found the photon distribution functions fm{n) for all modes. For the modes 
which do not exhibit squeezing, the PDF coincides with a time-dependent Planck's 
distribution, while the PDF in the distinct 'principal' modes differs from Planck's one, 
being more 'flat' for n ^ 1. 
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Appendix A. Some properties of the Bogoliubov coefficients p(^) 

The non-zero coefficients p^^^ for the 'principal' modes read 

ipn+p/2) ^ T{n + 3/2) ^n-m^m+n+l 

Ppm+p/2 r (m + 3/2) r (1 + n - m) 



13 



X 



F(n + 1/2, -m-1/2; 1 + 



n — m: K 



n > m 



(A.l) 



{pn+p/2) ^ (-l)'"-"r (m + 1/2) 
Ppm+p/2 



X 



r(n + l/2)r(l + m-n) 
F (^m + 1/2 , —n — 1/2 ; 1 + m — n ; , m > n 



(A.2) 



(pn+p/2) ^ (-l)"^r (m + 1/2) r (n + 3/2) «:«+™+iA^ 

P—pm—p/2 



Using the formula 



vrr (2 + n + m) 
X F (^n + 1/2 , m + 1/2 ; 2 + n + m; 

T{a + b + l) 



F{a,b]a + b+ 1; 1) 



r(a + l)r(6+l) 

one can find the asymptotics of the coefficients p^-* for k ^ 1 |T7 

Pttfir » 1) ^ sin[rr(m + j7p)] ^ ^ +„+2,/p^n-™ 

In particular, for p = 2 and for the odd ('principal') modes 

p'^2:^ir>>l)^^t}r_^a + ^^r-- 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



7r(2m+ 1) 

It is known that the hypergeometric function F{a, b; c; z) with 'half-integral' 
parameters a, b and an integral parameter c can be expressed in terms of the complete 
elliptic integrals 

rV2 da TT 



E(k) 



Vl — sin^ a 2 

7r/2 



..,11 

-F -, - ; 1; 
- ' 2 ' 2 



1 / ^ 2~ ^ r. f 1 1 

day 1 — K"' sm a = —F — , - ; 1; 

^ 2 V 2 ' 2 ' ' 



In particular, 



a)_2X{K) a) 2 

Pi Hn^ P-i 



(1) 



E(k) - k^K(K] 



P-3 



2A^ 

37rK 



1 - 2k'^] E(k) - k'^K(K) 



2 - k']E(k) - 2k^K(K 



(A.8) 
(A.9) 

(AlO) 
(All) 
(A12) 
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where k = yl^^fi? and A(k) = y/1 — 72^2 _|_ 

The general structure of the coefficients p^^^ (we confine ourselves to the case p = 2) 
is as follows 

(1) 2X^+H^) 

P2m+1 - 



TCK" 



P-2m-l 



fm E{k) + k^Qm K(/€) 

^ (k^) E(k) + K^Sm (/t^) K(/t) 



(A.13) 

(A.14) 



where fm{x), gm{x),rm{,x), Sm{,x) are the polynomials of the degree m which can be 
found from the recurrence relations (|2.6| ). 

The Taylor expansions of the complete elliptic integrals at k — (when n f» apr) 

read 



E(/t) 



64 



K(/€) 



TT 



, 1 2 9 
1 + -K^ H K 

4 64 



whereas their asymptotic behaviours at k ^ are given by the formulae 
K(fi;) fti In - + - In 1 H 

E(k) ^ 1 + - fin - - + ■ • • . 

In this case ~ a/ sinh(apr) and ln(l/«:) ^ apr. 



Appendix B. Calculation of integrals 



To calculate, for instance, the variance Ui we use equations (|3.7|) and (|A.1CI| ) and replace 
the derivative over r by the derivative with respect to k using the relation (if p = 2) 
dn = 2/3ft:^dr, where P = yr^^^pK^. We arrive at the equation 
dUi 2 



{[^2 (i - 2-f\'^) + 1 - 2(3 



-2k^{l - /5k)E(k)K(/€) + k'K^k) \ . (B.l) 
Let us consider first the case 7 = 0, when (3 = 1. Taking into account the differentiation 



rules 



dK(/t) E(/t) K(k) dE(K) E(k) - K(k) 



dK KK^ K dK K 

we may suppose that the factor k^ in the denominator of the right-hand side of equation 
( [B.l| ) comes from the derivative dK/dfi;. Thus it is natural to look for the solution in 
the form 

2 



-k'^k 



A{k)K\k) + B{k)K{k)E{k) + C{k)E\k) 



(B.3) 



where A{k), B{k) and C{k) are some polynomials of k. Putting the expression 
into equation ( [B.l| ) we obtain a set of coupled equations for the unknown functions 
A,B,C. Writing A{k) = oq + Ai{k), B{k) = bo + Bi{k), C{k) = cq + Ci{k) we 
determine the constant coefficients oq, &o and Cq by putting k = in that equations. 
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Then we obtain new equations for the functions Ai{k), Bi{k,) and Ci(k) and repeat 
the procedure. After a few steps we arrive at the equations which have obvious trivial 
solutions An = = C„ = 0. This confirms our hypothesis on the polynomial structure 
of the functions A{k), B{k) and C{k) and gives the final answer. The equations for 
the variances f/^, V^, etc. with > 3 can be integrated in the same manner, the only 
difference is that one should write instead of k in the denominator of the expression 
like ([B.3D . In the generic case 7 7^ we notice that the factor (3 can appear in the 
denominator of the expression ( [B.l| ) as a result of differentiating the function 
since d/?/dK = — 7^k//5. Therefore we split each function, A, 5,(7 in the '/5-even' and 
'/?-odd' parts like A = Ae{n) + i3{k)Ao{k). The equations for the 'even' and 'odd' 
coefficients turn out independent, and we solve them using the procedure described 
above. The equation (|3.21|) was integrated using the same scheme. 



Appendix C. Asymptotics of the minimal variance and purity factor 



For the initial diagonal density matrix (in the Fock basis), combining the equations 
( p.ll| )-( p.l3| ) and ( [4.5|) , we write the variances at r ^ 1 as (we omit the subscript fi) 



( U{t) \ 
V{r) 



I 2Fsin2(x/2) + / \ 



(ar + Z) 



J 



X = /i0. (C.l) 



2Fcos2(x/2) + (? 
\^ — F sin X + ^ 

The functions /, g and h are much smaller than F. At r ^ 00 these functions tend 
to finite limits which do not depend on the initial state, since they can be found by 
integrating equations ( p.ll| )-( [3rT^ ). Thus we have U + V = 2F + f + g, whereas 

{U - V)^ + = 4F^ + 4F[{g - /) cos x - 2/i sin x] + (/ - g)^ + 4h^. 

For F ^ f,g,h we have 



'{U - Vy + 4F2 = 2F + {g-f) cos x - 2/i sin x + 0{1/F) 

so the minimal variance u{t) ( ^.15| ) tends to the finite limit 

u{oo) = /cos^(x/2) +5(sin2(x/2) + /isinx 

which does not depend on Z, i.e. on the initial state. 

Analogously, UV — Y^ = 2Fu{oo) + 0{1) ~ r for r ^ 1. Consequently, the purity 
factor X asymptotically goes to zero as r"^/^. 



Appendix D. Photon distribution in the Gaussian state 

In the most compact form the information on the photon distribution f{n) in some 
mode (we suppress here the mode index) is contained in the generating function 

00 

G{z) = E /(^)^" 

71=0 
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For the most general Gaussian state it was given in ^ (for a single mode): 



where 



D 



zgi - z'g2 



90 



(D.l) 



l + zf + i (UV - Y^) (1 - zf + 2{U + 1/) (l - z^)] (D.2) 

D = l + 2{U + V)+4: (uV - Y^) = 4^(0) 
go = {P)\2U + 1) + {q)\2V + 1) - 4(p) {q)Y 

gi = 2{pf [u^ + Y^ + U + ^+ 2{qf (v^ + Y^ + V + ^ 
-A{P){q)Y{U + V + l) 

g, = 2{pf {u' + - i) + 2(^)' + - i) - ^^P)m{U + V) 

If {p) = (q) = 0, then the probability f{n) is expressed in terms of the Legendre 
polynomials: see equation (|5.1| ). In the generic case f{n) is related to the two- 
dimensional 'diagonal' Hermite polynomials ||38| : 

fin) = ^Hi^Hx,x*) (D.3) 



where 



J^o = f{0) = 2D-'/^exp{-go/D) 

V2{i2V - l){q) - 2Y{P) + z [(1 - 2U){p) + 2Y{q)]} 



X 



2{U + V)-A{UV -Y^) - 1 
and 2x2 symmetric matrix TZ has the elements 



■Rii = ^-(V-U- 2,y), Tli2 



1 - 4 (uV - Y' 



The two-dimensional Hermite polynomials are defined via the expansion 

1 



exp 



-aJZa + aJZx j = ^ 



m\n\ 



(D.4) 



' m,n=0 

where x = (xi,X2), a = (01,02). The properties of these polynomials were studied 
recently in |3^, In particular, they can be expressed as finite sums of the products 
of the usual (one-dimensional) Hermite polynomials. The corresponding formula for the 
probabilities reads 



f{n) = To 



E 



|-f^n-fc(Or 



(D.5) 



where 



A 



(f/ - \/)2 + 4^2, S = A[\JV -Y^]-\ 



Squeezing and photon distribution in a vibrating cavity 17 
{2V + l){q) - 2Y{p) + z [(1 + 2U){p) - 2Y{ 



[2D{V -U - 2iY)] 



1/2 



The photon distribution function ( ^.1|) can be simphfied in the long-time hmit 



r ^ 1, when the average number of created photons N" = n ^ {V + U)/2 exceeds 
1. Then the mean-square fluctuation of the photon number has the same order of 
magnitude as the mean photon number itself, ^/On ~ V^A/", and the most significant part 
of the spectrum corresponds to the values n 1. Using the Laplace-Heine asymptotical 
formula for the Legendre polynomial P7 



Pn{z)^ ^ ^ , n»l 



^n+l/2 



one can simplify ( |5.1|) for the fixed values of the invariant variances u and v as 

fin) - ^=^= f 7^) (D.6) 



-Kn(v — u) V2z; + 1 

provided the positive difference v — u is not too small. Another approximate formula 
can be used if f ^ 1 but m ~ 1: 



Squeezing and photon distribution in a vibrating cavity 



18 



References 

Dodonov V V, Klimov A B and Man'ko V I 1989 Phys. Lett. A 142 511 
Schwingcr J 1993 Proc. Natl. Acad. Sci. USA 90 958, 2105, 4505, 7285 
Barton G and Eberlein C 1993 Ann. Phys. (NY) 227 222 
Lambrecht A, Jaekel M-T and Reynaud S 1996 Phys. Rev. Lett. 77 615 
Dodonov V V and Klimov A B 1996 Phys. Rev. A 53 2664 
Johnston H and Sarkar S 1996 J. Phys. A: Math. Gen. 29 1741 
Mundarain D F and Maia Neto P A 1998 Phys. Rev. A 57 1379 
Dalvit D A R and Mazzitelli F D 1998 Phys. Rev. A 57 2113 
Schiitzhold R, Plunicn G and Soff G 1998 Phys. Rev. A 57 2311 
Janowicz M 1998 Phys. Rev. A 57 4784 
Ji J-Y, Jung H-H and Soh K-S 1998 Phys. Rev. A 57 4952 
Ji J-Y, Soh K-S, Cai R-G and Kim S P 1998 J. Phys. A: Math. Gen. 31 L457 
Lambrecht A, Jaekel M-T and Reynaud S 1998 Eur. Phys. J. D 3 95 
Golestanian R and Kardar M 1998 Phys. Rev. A 58 1713 
Dodonov V V 1998 Phys. Rev. A 58 4147 
Pendry J B 1998 J. Mod. Opt. 45 2389 
Dodonov V V 1998 J. Phys. A: Math. Gen. 31 9835 

Mendonga J P F, Maia Neto P A and Takakura F I 1999 Opt. Comm. 160 335 
Wu Y, Chan K W, Chu M-C and Leung P T 1999 Phys. Rev. A 59 1662 
Miri F and Golestanian R 1999 Phys. Rev. A 59 2291 
Wu Y, Chu M-C and Leung P T 1999 Phys. Rev. A 59 3032 
Dalvit D A R and MazziteUi F D 1999 Phys. Rev. A 59 3049 

Plunien G, Schiitzhold R and SofF G 1999 Los Alamos archive preprint |quant-ph/9906122 
Cooper J 1993 LEEE Trans. Antennas Prop. 41 1365 
Dittrich J, Duclos P and Gonzalez N 1998 Rev. Math. Phys. 10 925 
Dodonov V V, Klimov A B and Man'ko V I 1990 Phys. Lett. A 149 225 
Chizhov A V, Schrade G and Zubairy M S 1997 Phys. Lett. A 230 269 
Dodonov V V and Khmov A B 1992 Phys. Lett. A 167 309 
Dodonov V V, Klimov A B and Nikonov D E 1993 J. Math. Phys. 34 2742 
Moore G T 1970 J. Math. Phys. 11 2679 
Law C K 1994 Phys. Rev. A 49 433 

1995 Phys. Rev. A 51 2537 
Dodonov V V, Kurmyshev E V and Man'ko V I 1980 Phys. Lett. A 79 150 

Dodonov V V, Klimov A B and Man'ko V I 1993 Squeezed and Correlated States of Quantum 
Systems (Proc. Lebedev Physics Institute 205) ed M A Markov (Commack: Nova Science) p 61 
Dodonov V V, Man'ko V I and Polynkin P G 1994 Phys. Lett. A 188 232 

Dodonov V V and Man'ko V I 1987 Group Theory, Gravitation and Elementary Particle Physics 

{Proc. Lebedev Physics Institute 167) ed A A Komar (Commack: Nova Science) p 7 
Dodonov V V, Malkin I A and Man'ko V I 1975 J. Phys. A: Math. Gen. 8 L19 

1975 Int. J. Theor. Phys. 14 37 
Dodonov V V and Man'ko V I 1989 Invariants and the Evolution of N on- Stationary Quantum 

Systems {Proc. Lebedev Physics Institute 183) (Commack: Nova Science) p 263 
Dodonov V V, Man'ko O V and Man'ko V I 1994 Phys. Rev. A 49 2993; 50 813 
Dodonov V V and Man'ko V I 1994 J. Math. Phys. 35 4277 
Vourdas A and Weiner R M 1987 Phys. Rev. A 36 5866 
Schleich W and J.A.Wheeler J A 1987 J.Opt.Soc.Am. B 4 1715 
Dodonov V V, Klimov A B and Man'ko V I 1989 Phys. Lett. A 134 211 
Dodonov V V 1994 J. Phys. A: Math. Gen. 27 6191 

Gradshtein I S and Ryzhik I M 1994 Tables of Integrals, Series and Products (New York: Academic) 



Squeezing and photon distribution in a vibrating cavity 19 

[45] Prudnikov A P, Brychkov Yu A and Marichev O I 1986 Integrals and Series. Additional Chapters 
(Moscow: Nauka) 

[46] Erdelyi A (ed) 1953 Bateman Manuscript Project, Higher Transcendental Functions (New York: 
McGraw-Hill) 

[47] Szego G 1959 Orthogonal Polynomials (New York: American Mathematical Society) 



Squeezing and photon distribution in a vibrating cavity 
Figure captions 



20 



Figure 1. The minimal variance wi, the maximal variance di, the purity factor the 
mean photon number Mi and Mandel's parameter Q of the fundamental mode jjl = 1 
versus the universal parameter k (2.8). 



Figure 2. The photon distribution function in the fundamental field mode /i = 1 
for r = 5 and 7 = 0, p = 2 (points connected with solid lines). Points connected 
with dashed lines correspond to the Planck distribution with the same mean photon 
number. 



